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Abstract. In this paper we study the construction of a discrete solution for a hyperbolic system of partial 
differentials of the strongly coupled type. In its construction, the discrete separation of matricial variable 
method was followed. Two separate equations in differences were obtained: a singular matricial and the 
other one a Sturm Liouville vectorial problem, which by the superposition principle yield a stable discrete 
solution. 
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Resumen. En este trabajo se construye una solucion discreta a un sistema en derivadas parciales hiper- 
bolico de tipo singular fuertemente acoplado. En su construction se siguio el metodo de separation de 
variables matricial discreto, obteniendose dos ecuaciones en diferencias por separado: una matricial sin- 
gular y la otra un problema de Sturm Liouville vectorial, las cuales mediante el principio de superposition, 
producen una solucion discreta estable del sistema. 

Palabras claves: Sistema hiperbolico, contorno discretizado, problema mixto y de Sturm Liouville, solucion 
estable. 

1. INTRODUCCION 

Los sistemas hiperbolicos en derivadas parciales con coeficientes matriciales aparecen a menudo en 
la modelacion de situaciones tales como el estudio de los de procesos de calefaccion en microondas |6, 
[T3l , optica , cardiologia [T7] , fenomenos sismicos en medios elasticos H5] , entre otros. En este articulo 
consideraremos sistemas hiperbolicos en derivadas parciales con coeficientes matriciales constantes, 
definidos en el rectangulo [0, 1] x [0, T], los cuales caracterizan la ecuacion de onda homogenea: 

(1) Eu tt (x,t)-Au xx {x,t)=0, xe [0,1], 0< t < T, 

A 1 u{0,t) + A 2 u x (0,t)=0,0<t<T, 
lJ Bxu{l, t) + B 2 u x {l, t) =0< t < T, 



m u{x,0) = f{x), xe [0,1], 

l J u t {x,0) = g(x), xe[0,l], 

donde u = (ui,...,u m ) T , u tt = {u ut ,...,u m tt) T , u xx = {ui xx ,...,u mxx ) T ; f = (fi,...,f m ) T ,g = (gi,...,g m ) T e 
C m , A, A\, A 2 , B\, B2, E e C mxm , en su estudio d2) es llamada condicion de contorno fuertemente aco- 
plada y (13) el problema mixto asociado a CD - ©. El sistema QJ - (3) es tambien llamado singular si E 
es singular de lo contraio es conciderado no singular. Una gama de casos no singulares de CD - © han 
sido ampliamente tratados (vease HU, ED, E], QT|, El)- Si bien la literatura sobre su caso singular es 
muy poca \W\ estudia una mejora a [16] la cual soluciona un caso particular de CD _ 0. Este articu- 
lo se propone formalizar y generalizar los resultados de |10 | analizando la flexibilidad de las hipotesis 



Date: 19 de enero de 2013. 

2000 Mathematics Subject Classification. 35A09, 35L05. 



2 



MANUEL J. SALAZAR, EDISON E. VILLA 



que se asumieron para la obtencion de una solution al sistema de ecuaciones mencionado. El orden a 
desarrollar en el presente articulo es el siguiente: 

-Section 2: Se presentan definiciones y teoremas generates de las teoria de sistemas matriciales sin- 
gulares y de las ecuaciones en diferencias finitas, en particular demostraremos 3 proposiciones indis- 
pensables para la solution de la version discreta de CQ - © . 

-Section 3: Luego de discretizar el sistemaCD - ©, aplicando el metodo de separation de variables 
matricial discreto, se construye una solution a la ecuacion de onda considerando su condition de con- 
torno y se enuncia el teorema que resume el resultado obtenido. 

-Section 4: Se estudia el problema mixto y la estabilidad de la solution general. Con base en el teorema 
anterior enunciamos el segundo que garantiza la existencia de una solution discreta y estable de (D - @ . 

-Finalmente complementaremos lo obtenido con un ejemplo ilustrativo y enunciaremos algunas ob- 
servaciones y conclusion finales. 

2. PRELIMINARES 

En esta section damos algunos resultados y notaciones que seran de utilidad para las siguientes 

Definition 1. SeaAe c mxm i una matriz, la cual denotamos por A G , es llamada inversa generalizada de 
A si veriflca la condition 

AA G A = A 

Teorema 2. {Teorema de Mitra) Sea A e c mxm y a g una inversa generalizada para A. Sea b e C m , 
entonces 

Ax = b tiene solution o AA G b = b 
y cualquiera de sus soluciones tiene la forma 

(4) x = A G b + [l-A G A)z, zeC m . 

Demostracidn: Vease H2] . 

Definition 3. Sea W un subespacio de C m y A una matriz en c mxm , se dice que W es invariante por la 
matriz A, si veriflca que AW c W {i.e Vx e W => Ax e W) . 

Teorema 4. Sean A, B, B G e c mxm i entonces B A [I - B G B) = siiKer [Bi) es un subespacio invariante de 
A. 

Demostracidn: Vease H2] . 

Definition 5. (La Inversa Drazin) Sea A £ c mxm tal que ind(A) = k, la inversa Drazin de A (A D ) es la 
unica matriz que cumple 

1) A D AA D = A D 

2) A D A = AA D 

3) A k+l A D = A k 

Teorema 6. Este teorema consta de 3 numerates: 
1. Si A tiene la descomposicion canonica de Jordan 

'CO 



(5) A = T 



N 



donde C e C pxp y T e c mxm son matrices invertibles, N e C qxq es una matriz nilpotente con indice k con 
[p, q) en N 2 cumpliendo p + q = n, entonces la inversa Drazin de A, A D , se obtiene de la siguiente forma 



(6) A D = T 



C" 1 ] - 







T 
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2. Para cierto polinomio Q (x) , La inversa Drazin de A puede ser escrita como 

A D = Q (A) 



A G =T 



1 C" 1 ) r _! 



3. Si en la descomposicion canonica de A en © N = 0, entonces 

r 1 o 



Demostracion: Vease |2] y [12] . 

Definicion 7. Sea A £ £ mxm i f:C—*C, decimos que f pertece al espacio 9 {A) si existe un entorno V del 
espectro de A, a {A) sobre el cualfes anaMtica. 

Teorema8. [Aplicacion espectral) Sea f e 9(A), entonces 

f{a{A))=a{f{A)), 

donde f {a (A)) = {f (A) : A e a (A) } . 
Demostracion: Vease (3 . 

m 

Definicion 9. Seax £ C m . deflnimos la norma vectorial \\ \\ \ por: \\x\\ i = || (x\, x m ) \\i= £ |x;| , yla norma 

i=l 

matricial por ella inducida como 

(7) \\A\\ X =sup{||Ax|| 1 :xeC m > ||x||i = l}, paraAeC mxm 

si las columnas de A son a\, a.2,—, a n para las cuales deflnimos sus componentes aij, i = I,..., n, para 
j = 1, n respectivamente, entonces cohincide con 



|| AH i =max|^ \atj\ : j = 



Teorema 10. Sea A e c mxm , / e9(A), lafuncion matricial f (A) se define por la matrix. 

f(A) = P(A) 

donde P(x) es un polinomio conocido de grado menor que el del polinomio minimal deA. 
Demostracion: Vease . 

Teorema 11. Sea A e £ mxm talque a — max {\z\; z e a (A)} , para e > existe una norma matricial en 
C mxm talque \\A\\<a + e. 

Dada la variedad de tipos de estabilidad manejadas en la literatura de los esquemas de diferencias fi- 
nitos numericos asociados a ecuaciones en derivadas parciales, nosotros hemos optado por la siguiente 
definicion: 

Definicion 12. \T%\{Estabilidad). Consideremos fijo el par ordenado {h, T) £ R + x |R + . Si se tiene una 
malla de puntos en el piano (x, t) con paso espacial h y paso temporal k > talque M es un numero 
natural. Si U es unafuncion definida en dicha malla y U[ih, jk) representa su valor en el nodo [i, j) , se 
dice que U es estable, cuando U[ih, jk) permanece acotada independiente de los valores dekyM tal que 
kM=T 

Definicion 13. (Problema de Sturm Liouville Discreto) El problema de valor en la frontera formado por 
la ecuacion en diferencias 

(8) k[p(i-l))Au(i-l) + q(i)u(i)+Ar(i)u(i)=0, 0<i<N 
con condiciones de frontera 

(9) u{0) = au(l), u{N) = Pu(N-l), 
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en las que lasfunciones son escalares se llama problema de Sturm- Liouville discreto (PSD). 
Donde A es el operador de diferencia deflnido por 

Ap (z -l) = p (i) -p{i-l), A 2 p (i) = A (Ap (i)) . 

en la ecuacion en diferencias A esun pardmetro, ay (5 son constantes conocidas, lafuncion p (z) esta 
deflnida enO < z < N- 1, p [i] , q{i) y r (z) enO < i < N; ademds p (z) > yr [i] > 0. 
Al efectuar la operacion A en la ecuacion en diferencia^, encontramos para < i < N 

p (z) u {i + 1) - [p (z) + p{i- 1)) u (z) + [q (z) + Ar (z)) u (z) + p (z - 1) u (z - 1) = 

denotando s(i) = p(i) + p(i -1) - q(i) ,0 < i < N, se tiene 

(10) -p{i- l)u[i - 1) + s(z) iz(z') -p{i)u{i + 1) = Ar (z) u(i) , < i < N 
en esta ultima ecuacion para los dos casos k = 1, N- 1 tenemos las ecuaciones 

-p(0)M(0) + s(l)M(l)-p(l)M(2) = Ar(l)zz(l), 
-p(N-2)u(N-2) + s(N-l)uiN-l)-p(N-l)uiN) = Ar (N -1) u(N -1) , 
las cuales teniendo en cuenta @ reescriben como 

(11) s(l)zz(l) -/?(!)" (2) = Ar(l)u(l), 

(12) -p(N-2)u(N-2) + s(N-l)u(N-l) = Xr{N-l)u(N-l), 
dondesa) = sa)-ap(.0)ys(N-l) = s{N-l)-Pp{N-l). 

Con las ecuaciones {77J para z = 0, Q2 para z = iV— 1 y 03 para 2<i< N-2 se producen N- 1 ecuaciones 
equivalentes a @ - @ , /as cuales describen un problema generalizado de valores propios matricial con la 
forma siguiente 

a(l) 
u{2) 



u{N-2) 
[ u(N-l) 

u(l) 
u{2) 

u{N-2) 
u(N-l) 

Ya que p (z) > 0, < i < N— 1; r (z) > 0, < z < N de acuerdo a la teoria espectral para matrices simetricas 
nosotros asumimos los siguientes teoremas. 

Teorema 14. El PSD tiene exactamente N -I eigenvalores A/, < I < N; a cada uno de los cuales les 
corresponde las eigenfunciones vi{i), < i < N deflnidas para < i < N. 
Demostracion: Vease Q] 6 HU . 

Teoremal5. Sean Xi y vi{i) losN-l eigenvalores y egenfuciones del PSD, entonces el conjunto 

{ Vl (i):l = l,...,N-l} 
es ortogonal con respecto al producto interno 

N-l 

£ r (z) (z) v v (z) =0, Vp ^ v. 
i=l 



(13) 



s(l) 

-pa) 






-pa) 

5(2) 
-p(2) 





-P(2) 



r(l) 









r(2) 




s{N-3) 
-p{N-2) 



-p(N- 
s(N-1 






r{N-2) 
r(AT-l) 
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Demostracion: Vease Q] • 

Teorema 16. Sea {vi(i) : < / < N} el conjunto de eigenfunciones en el teorema anterior y u (z) cualquier 
funcion deflnida para < i < N , entonces u (z) puede expresarse como una combination lineal de dicho 
conjunto, es decir: 

N-l 

u (z) =Y,ciVi (z) , < i < N 



donde cada c\ puede ser determinada por 

LfS^rii) vi(i)u(i) 

ci = 



0<1<N 



L-^rm v\{i) 
Demostracion: Vease HI . 

Las siguientes 3 proposiciones las enunciamos con una notation pertinente para la demostracion de nues- 
tro teorema fundamental (l22t . 



Proposicion 17. \Sean ay/3 escalares conocidos, N e N y X un eigenvalor del PSD 

(14) h(i + l)-(2-X)h(i) + h(i-l) = 0,0<i<N, 

(15) h(.0) + aN(h(l)-h(.0)) = 0, 

(16) h{N) + pN{h{N)-h(N-l)) = 0, 

con eigenfuncion v{i). Entonces v (z) es solucion del PSD Q3l - GUl si y solo si 

H{i) = v(i)R 

con R vector arbitrario en C m , es una solucion del problema de Sturm Liouville discreto vectorial (PSDV) 

(17) H{i + l)-(2-X)H(i) + H(i-l) = 0,0<i<N, 

(18) H(0) + aiV(H(l)-H(0)) 

(19) H(N) + pN{H{N)-H(N-l)) 

en el que H{i) es una funcion vectorial discreta enC m . 
Demostracion: 

Supongamos que H (i) satisfaceelPSDVn7k-[Wb donde H(i) = [ h\{i) h.2{i) ... h m {i) ) siendoh^ii) 
funciones escalares discretas definidas para < i < N. Luego en virtud de la igualdad entre vectores tene- 
mos que cada una de sus componentes (z) A; = 1, m satisfacen el PSD fZ4t - (76| . Entonces al hallar 
una solucion para una componente, ella lo sera tambien para las demds. Sin perdida de generalidad de- 
notemos ht {i)-h (z) . Si en O deflnimos: 

q{i)=0, r(z') = 1, p(i) = I, 

de G3 - (O se obtiene 

2 ~ aN -10 

1-aN 1 u u 

-12-1 : 



0, 
0, 



(20) 







: 2 -1 

-1 ^ 

Por el teorema O , (20} tiene losN - \ eigenvalores 







' MD ' 


h{2) 




h{2) 




= X 




h(N-2) 




h(N-2) 


h(N-l) 




h(N-l) 



(21) 



{Xi 



N-l 
1=1 



6 MANUEL J. SALAZAR, EDISON E. VILLA 

con eigenfunciones asociadas 

(22) {Viii)}?^ 1 , < I < N, < i < N, 

vi = {vi(l), v l {2),„Vi(N-l)) T , 0<1<N. 
Luego, las soluciones del PSD QUi - O tienen la forma 

vi(i) = V\ (Z) D 2 , con D2 arbitrario en C e i = 1, ...N - 1. 
Como Zo anterior se veriflca para cada componente de H{i) concluimos que 

(23) H(i) = (hi (i) , fe w (Z)) = y z (i) i?, i? arbitrario en C m 
es solution del PSDVfLTl - 03 • 

Aplicando el teorema fZ4t cZaro aflrmar que el PSDV(T7l - O fZene exactamente N-l soluciones de 
la forma Q?3t donde v\ (i) son las N-l eigenfunciones de PSD fZ4t - (TT6b determinadas en f22li . 

Proposition 18. : SeaH(i) funcion vectorial discreta en C m , que veriflca 

(24) H(Z + l)-(2-A)H(Z) + H(Z-l) = 0, 0<i<N, 

(25) Ai#(0) + A 2 iV(H(l)-#(0)) = 0, 

(26) BiH{N) + B 2 N(H(N) - H(N- 1)) = 0, 

con Ai, A 2) Bi ) B 2 eC m><m y 

(27 > g («-«=(m-S)- 

talque rang[G[a, /3)) <myRe KerG[a,f>). Entonces para cada A/, < Z < iV eigenvalor del PSD fZ4t - 

H/(i) = i^(Z)i?, 0<i<N 

son N-l soluciones de f2it - (|25t . 
Demostracidn: Sea 

Hi(i) = vi (Z) R0 < i < N, Re KerG{p,p) 
con H{i) asi definida, f24t - (|25t es equivalente a 

(28) (i/(i + l)-(2-A)i/(i) + i/(i-l))i? = 0,0<i<N, 

(29) (Aii/(0) + A 2 iV(i/(l)-i/(0)))i? = 0, 

(30) (Bii/(N) + fl2N(i/(iV)-i>UV-l)))i? = 0, 

como i; (Z) veriflca QHl - O . ft'erce gue |[28) - J30i , es equivalente a 

[v (Z + 1) - (2- A) v (i) + v(i-l))R = 0,0<i<N, 

(aAi-A 2 )(v(l)-v(0))R = 0, 
[PB 1 -B 2 )(v(N)-v(N-l))R = 0. 

Como ademds ReKerG[a,p), se concluye que Hi(i) = v\ (Z) R, satisface (24& - (251 ■ 

El teorema y proposition siguientes los enunciamos con una notation pertinente para la demostracidn de 

nuestro teorema fundamental. 

Teorema 19. Sean I, Ae c mxm donde I es la identidad, yG{j) funcion vectorial definida para i = 0, ...N; 
supongamos ademds p escalary consideremos la siguiente ecuacion en diferencias matricial 



G [j + 1) - (2/ + pA) G [j] + G[j-1) = 0, j>0.G [j) e C 



mxm 
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Entonces su solution general es 

G[j) = [P + (A)] j P 1 + [P-(A)] j Q l , j>0, h,l 2 eC m , 

donde P+ (x) y P- (x) son los polinomios de grado mas pequeho que el grado del polinomio minimal de A 
tal que 

P + (A) = I+^A +] J(i+^a) 2 -I,P- {A) = I+^A-^I+^A) 2 -I. 
Demostracion: Vease [3|, o tambien |16|. 

Proposition 20. Con las mismas hipotesis del teorema anterior, consideremos la siguiente ecuacion en 
diferencias matricial 

(31) EG{i + \)-{2E + pA)G{i) + EG{i-\) = 0, ; > 0. 

donde EeC mxm es singular, supongamos ademds que existe y eC tal que (y£ + A) es no singular, deflni- 
namosE = {jE + A)~ 1 E, A = {XE + A)~ l A. Entonces su solution general vienedada por 

G{j) = Z J EE D h + Z Y j EE D l 2 , con h, l 2 arbitrarios en C m . 
donde Z = [P+ {E D A)} EE d , Z x = [P- [E D A)] EE d 

con 

P+ (E D A) = 1 + ^E D A + )J[l+^E D A} 2 -I, P- (E D A) = 1+ ^E D A - J\l+-E D A^ -I 

Demostracion: Dado que existe y con {jE + A) invertible, sean E = {jE + A) ~ 1 E, A = [XE + A) ~ 1 A, luego la 
ecuacion OTT i es equivalente a 

(32) EG{j + l)-(.2E-pA)G[j) + EG[j-l) = 0, j>0. 
Considerando la descomposicion canonica de E, obtenemos: 

en lasqueCeC pxp es matrix invertible, N e C qxq esuna matriz nilpotente de indice k con p, q cumpliendo 
p + q = n. Luego haciendo 

{ Hi) 



[d(j) }' 



(34) TG{j) 
(T3TT) es equivalente a los dos siguientes sistemas: 

(35) Ch[j + l)-[2C + p[l-rC)]h[j) + Ch{j -l) = 0, 

(36) Nd(j + l)- [2N + p{l-rN)]d[j) + Nd{j -1) = 0. 
Como C es invertible se obtiene que (33 se veriflca si y solo si 

h [j + 1) - [21 + pC" 1 (J -TC)]h [j) + h [j - 1) = 0, 

la cual por el teoremalHh tiene la solution general 

h [j] = [P + (D)V Pi + [P. {D)]J Qx, ; > 0, h, l 2 e C k . 

donde D = C" 1 (l - yC) y P+ (x) y P- (x) son los polinomios de grado mas pequeho que el grado del poli- 
nomio minimal de D tal que 



p + {D) = i + ^d + J(i + ^d) 2 -i,p.{D) = i + ^d-J(i + ^d) 2 -i 
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Como para N suflciente mente grande 2N + p[l- yN) es invertible ykesel indice de nilpotencia de N, 
por induction puede demostrarse que la ecuacion OB tiene la solution trivial d [j] = 0, j > 0, luego con- 
siderando 

(37, e=t-'[ c "jr.^r-^ 1 °)r 

y haciendo uso de las ecuaciones OH , 03 y OH , se concluye que 03 tiene la solution general 

(38) G(j) = Z j EE D h + Z 1 J EE D l 2 , con h, l 2 arbitrarios en C m . 

(39) dondeZ = [P+ [E D A)) EE D , Z x = [P- [E D A]] EE D 

3. SOLUCION DISCRETIZADA DE LA ECUACION DE ONDA 

Retomemos nuestro problema central 

(40) Eu tt {x,t)-Au xx {x,f) = 0, xe[0,l], t>0, 

(41) Aiu(0,t)+A 2 u x (0,t) = 0,t>0, 

(42) Biu{\,t)+B 2 u x (Y,t) = 0,t>0, 

(43) u(x,0) = f(x), xe [0,1], 

(44) u t (x,0) = g(x), xe [0,1], 

En el que asumimos ademas la hipotesis: que existe y g C tal que (jE + A) es no singular, y asi definimos 
E = (y E + A) ~ 1 E , A = ( XE + A) ~ 1 A con la condicion: 

(45) z ^ para algun zea [E D A) c C 

En el esquema de diferencias correspondiente a |[40l - |[44ll dividimos el dominio [0, 1] x ]0, oo[ en rectan- 
gulos de lados Ax = h, At - k, luego al introducir coordenadas de un punto tipico de la malla [ih,jk) , 
representamos el valor u[ih,jk) por U[i,j). Si aproximamos las derivadas parciales de u utilizando 
diferencias avanzadas para las derivadas primeras y diferencias centrales para las derivadas segundas. 
Ellas tienen la forma 



(46) u t (ih,jk) 

(47) u x (ih,jk) 

(48) u xx {ih,jk) 

(49) u tt (ih,jk) 



U{i,j + l)-U{i,j) 
k 

U{i + l,j)-U{i,j) 
h 

U{i + l,j)-2U{i,j) + U{i-l,j) 
h 2 

U{i,j + l)-2U{i,j) + U{i,j-l) 

k 2 



Al sustituir |[46l - |[49l en J40l l - (l44l l , para N un entero positivo con h = ll Ny r = k/hO < i < N;se obtiene 
la siguiente representacion aproximada de 0OJ - 04} 



(50) 



r 2 A{U{i + \,i)-U{i-\,i))+2{l-r 2 A)U{i,j)- 
{U{i,j + l) + U{i,j-l)) = 0, 0<i<N, j>0 
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(51) A 1 U(0,j) + NA 2 [ua,j)-mO,j)) = 0, j>0, 

(52) B l U{N,J)+NB 2 [mN,j)-mN-l,j))=0, j>0, 

(53) U(i,0) =F(i), 0< i< N, 

ua,i)-mm n . 

(54) = G{i), 0<i<n. 

k 

Busquemos soluciones no trMales a la ecuacion de onda discretizada ([50]) . Para esto usamos el metodo 
de separation de variables discreto matricial suponiendo U [i, j) de la forma 

(55) U[i,j) = G[j)H(i), H(i) EC m . 

Para r = klh, imponiendo que {U{i,j)} verifique (|50), resulta 

r 2 AG{j)[H(i + l)-H{i-l,j)}+2{E-Ar 2 )G{j)H{i)- 
1 J E[G(j + l) + U[j - l)H(i)] = 0, 0< i<N, j>0 

Tomando peU arbitrario y puesto que 2 [E - Ar 2 ) = - r 2 A \l + ^ j + (2£ + p A) , ([5GJ toma la forma 

(57) r 2 AG[j)[H{i + l)-A(2 + ^]H(i) + H(i-l) - 
[EG[j + l)-{2E + pAG[j)+EG[j-l)]H{i)=0, < i < N, j > 0. 

La ecuacion J57J se verificasi {H(/)} y |G(;)} satisfacen las ecuaciones en diferencias 

(58) H(i + l)-(2 + -^H(i) + H(i-l)=0, 0< i< N, 

(59) EG(j' + l)-(2£ + pA)G(j)+£G(7-l)=0, ; > 0. 

Si en ([58) definimos p = -r 2 Xi para cada A; eigenvalor del PSD {14} - CGS , de la proposition {T7J , sabe- 
mos que si v\ (i) es el eigenvector asociado a A/. Entonces: Hi (i) definida como 

(60) Hi(i) = v t {i)R, ReC m -{0}, 0<i<N, j>0. 
satisface ll58l . 

Por otra parte usando la propocision GUJ, para p = -r 2 Xi, la ecuacion |[59), tiene solution dada por 

(61) G{j) = Z ] E^h + Z^ EE D l 2 , con h, / 2 arbitrario s en C m . 

(62) donde Z = [P+ {E D A}] EE D , Z x = [P. [E D A]] EE D 
asi, una solucion de l[50| es de la forma 

(63) U{i, j) = Z j H u (Z) + Z j H 2l (i) , < i < N, j > 0. 
donde 

(64) H u ii) =P t v ii) , H 2l (*) = QiV (i) con P h Q l e C m - {0} . 

3.1. Condition de contorno discretizada. Las soluciones U[i, j) de lEOj encontradas en ([63), se escri- 
ben de acuerdo a lEU , como sigue 

(65) Ui [i, j) = Z Q i P l vi (/) + Z^Qi vi ii) , < i < N, j > 0. 

donde v t (i) es alguna eigenfuncion del PSD {Hi - {H y P h Qi e (C m - {0}) . Sea G (a, /3) definida en {27) , 
por la proposition {18) , dejando ; fijo, {65| satisface (ED - {52J si P/ y Q/ en los sistemas de ecuaciones 

Z j (Aiv(0) + A 2 N(va)-vm)Pi = o, 
Z j (B l v(N) + B 2 N{v VSl)-v(N-l)))Pi = 0, 
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y 

Z x j (A 1 v(0) + A 2 N{va)-v(0)))Qi = 0, 

Z^iBivm + BzNivm-viN-mQi = 0, 

pertenecen a KerG(p, ft) n (C m - {0}). Con lo anterior Ui [i, j) en l[65) es solution de l[50j 
Los resultados obtenidos hasta el momento, los resumimos en el siguiente teorema. 

Teorema 21. Supongamos fix), g(x) funciones en C m , A, A\, A 2 , B\, B 2 ; E e c mxm , E matriz singular 
con y tal que (yF + A) es invertible, A, E, E D e C mxm deflnidas por lO, C33 y que veriflcan la condi- 
tion 03; sean ademds (a,/3) e U 2 que veriflcan lO - G3. G(a,/3) representada por (23 que cumple 
rang [G [a, ft)) < m. Entonces una solution de - lO esta dada por: 

N-x . 

(66) U{i,j) = £ [Z J EE D P l + Z{EE D Qi) v x (z) ,0<i<N,j>0 

donde Z y Z\ estdn deflnidas por (62| , l\ , l 2 e KerG [p,, /?) n (C m - {0}) y v\ (z) son lasN-l eigenfunciones 
asociadas al PSD O - (O determinadas segun la proposition 03 • 

4. Problema mixto 

Analizamos ahora las condiciones a imponer para que U{i,j) en (l66l l verifique el problema mixto 
(153) - . Para esto al sustituir J66J en dichas expresiones, se deben cumplir las siguientes dos ecuacio- 
nes: 

IV- 1 

(67) F(z) = £ ( ^-P/ + ^Q/) (0 



(68) G(i) = 



Z=l 

Lg" 1 [Z EE D P l + Z x EE D Qi) vi (z) - F (i) 



jfc 

con lo anterior, segun dTHJ , cada componente de F (z) y G (z) tienen su representation discreta respecto 
a {f / (z')}^\ luego las ecuaciones J67} y d68j son equivalentes a 



(69) EE u Pi + EE u Qi 



L i= i "zto 

(70) Z E£ u P l + Z 1 EE u Q l 



fy++ Dn * *f, D „ If-"/ ikG{i)+F{i)) 



que es un sistema matricial cuyas incognitas son los vectores Pi y Q/ los cuales son soluciones de las 



ecuaciones 



r 7 n rz 7 1&^p lf = - 1 1 ^q)[fcG(z) + (Z 1 -/)F(z)] 



(72) (Zi-^)M D Q 



I^ 1 vdiY 
.. /V) If," 1 ^(z)[(Z -/)F(z)-fcG(z)] 

l£>z(*) 2 



Sin perdida de generalidad, suponiendo que a (D) no contiene a 0, por el teorema (|20j , se tiene que 
a{D)=a (E D A) - {0} , ademas 

(73) a (P + (D)) = a (Z ) - {0} , a (P + (D)) = a [Z x ) - {0} 



SOLUCIONES DISCRETAS PARA SISTEMAS MATRICIALES EN DERIVADAS PARCIALES. 
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o sea 



luego 



a{P ± {D)) = 1 1 + 1<2 + ^[1+ |rf) 2 -l:dea (EA D ) - {0} j 



o- (P_ (D) - P + (D)) = |2^/(/ + |d) 2 -l :deff | 

y asi P_ (D) - P+ (D) es invertible si se verifica 

(74) pd(l + ^d)^0 

Luego tomando en (74) p suficientemente pequeno, utilizando la inversa de P_ (D) - P + (D) se pueden 
emplear los numerales 1 y 3 del teorema (2) , para hallar ((Zi - Z Q ) EE D ) G , la cual verifica: 

((Z! - Z ) £P D ) {{Z X - Z ) EE D f = EE D 

Asi, usando el numeral 2 del teorema (2) , se sigue que EE D conmuta con Z y Z : y por lo tanto los 
sistemas ED y £72) son consistentes si cumplen la condition 

(75) EE D F{i) = F(i), EE D G{i) = G(i) 

teniendo presente lo anterior, por el teorema (2), las ecuaciones GD y (72) admiten las siguientes solu- 
ciones 

(76) P/ = [(Z 1 -Z )EE j , 

2-i=l y /M 

(77) Q/ = [(Z 1 -Z )££ j . 

Con las condiciones 

(78) {G(i),F(i), 0<i<N}^KerG{n,P), 

G{n,P)E D A(l-G{Li,P) + Gln,P))=0, 

por el teorema (4), KerG{fi,/3) es un subespacio invariante de (£ D A). Por lo tanto con (75) P/ y Q/ 
cumplen 

EE D Pi = P h EE D Qi = Qi, {P h Q u O<i<N}czKerG{fi,p). 

lo que garantiza la consistencia de la condition de contorno discretizada (TSTb - (52l junto con el proble- 
ma mixto (53) - (54) • Luego, una solution de (50) - (54) viene dada por 



N-l. . . 
(79) U{i,j) = £ [Z J EE D Pi + Z{EE D Qi) v x (z) , 



N-l 

E 

i=i 

con Pi y Qi determinados por (76) y (77) . 
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4.1. Estabilidad. Para el analisis de la estabilidad de la solution C79), suponemos que F [j] y G [j] son 
acotadas. Por los teoremas C[0) y ©, de d73), el espectro de las matrices Z Q y Z Y en d) comprende los 
valores 



dado que 



J+fd + \/(/ + £d) 2 - J 



1, e a [EA D ], considerando que p puede verificar 
p < max{\Xi\; \\, 1 < I < AT- 1} x r 2 , 



(80) 

por teorema CCD se obtiene 

(81) ||Z || < ||Z 1 ||<l + 0(fc) ) ||Zo-/||<||Zi-/||<0(fc) 

|(Z -Zi)M D || < 0{k~ x ), jfc^O. 

para alguna norma matricial || || en c mxm . Usando ([TP , f72l . ED y el acotamiento de F(i) y G(z), con- 
cluimos que 



(82) 



|P/|| = 0(1), ||Q/|| = 0(1). 



Z ] 



permane- 



De lo anterior, se sigue que {U(i, j)} en C79) permanece acotada, si los numeros 

cen acotados cuando j — • oo, A; ->■ 0, 1 < j < M, Mk = T. Notese que de || Z 1| < 1 + O {k) se obtiene que 
para alguna constante positiva S, || Z 1| < 1 + kS, entonces, para 1 < j < M se calcula 



(83) 



Z ] 



M ^ „MO(Jt) ^ n MkS = e TS 



\Z \\ J <{l + (k)) J < (1 + O (fc)) M < e MUW < e 



. Asi, considerando |[82|, ([83j y L = max{|| 2//|| : 1 < I < N- 1], se sigue que la 



lo mismo ocurre para 

solution definida por l(79j , CHI y £77| es estable, es decir: 

I U [i, j) || = O (1) , fc — ► 0, h = — paso temporal fijo 
1 < z < N - 1 j — oo, con Mk = T 
En resumen, el siguiente resultado ha sido establecido. 

Teorema 22. Con las mismas consideraciones del teorema anterior, seanG(i) yF{i) acotadas quesatisfa- 
cen (73 y (IB con p veriflcando (EH y (8^ ; entonces una solucion estable del problema discreto (H^i - lO 
esta dada, por 



(84) 



N-l . 

U[i,j)= £ (Z^E^Pt + Z^E^QAvdi) 
1=1 



dondePi, Qi son como lo indican las ecuaciones y tZ3 
4.2. Ejemplo. Consideremos las matrices en C 3x 3 

e 8 
7 





f 





o 1 




, A = 





8 





, {e, j, a} eC-{0} 


t 


I o 





o- ) 





a 2 \ a 22 b 2 

«31 «32 b 3 j 



,A 2 



pan pa u ci 
\ia 2 \ pa 22 c 2 
\ pa 3l pa 32 c 3 j 



, B 2 = r)B x 
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13 





( /i(0 ) 




' gi m * 


F(i) = 




, G(i) = 






{ o ) 




{ o j 



escogiendo A ^ tal que Ay + 8 ^ 0, se calcula 

{XE + A)~ l 



( 1/A 8/Xe \ 
yl{Xy + 8) 



( \IXe -8lea " 

l/(Ay + 5) 

1/cr , 

/■ 



81 (Ay + 5) 



-o7e \ 








\ o 









o ) 


I o 







1 


J 










f A 





8Xle ' 




' 1 8/e ) 






f 












D _ 


o (a t 


+ 8)1 j 





, EE D = 


1 


, E D A = 





8/j 











I o 





o ) 




k ) 






I o 








J 


si en Q~4) - Q6} escogemos, a 




= t], entonces 


















f ° 








o 1 


T 








f ° 




























, G{a,p) + = 


\l{dl + dl + 


dl) 





















d2 d-i 





o ) 










k di 


d 2 


d 3 









donde di = \ib{ - q, i = 1, 2,3 y se comprueba ademas que: 

o-(£ D i) = {0, 1, 3} , ran (G (2,1)) <3, G(a,/3) F(i) = 0, G (a, jB) G(i) =0, 
EE D F{i)=F{i), EE D G{i) = G{i), G[iu, /3) E D All - G[[J,, + G(fi, fi)^ = 0. 

Luego, si p es suficientemente pequeno, las hipotesis del teorema (22j pueden ser verificadas y por lo 
tanto una solucion de (50) - d54) es de la forma 



N-i . 
U{i,j)= £ [Z l >EE D P l + ZiEE D Q l )vi(i) 
l=i 



deonde Pi y Q/ son como lo indican las ecuaciones J76j y J77J y v\ (i) son las eigenfunciones descritas 
en la proposicion dTTJi . 

4.3. Observacion. 

1 . El procedimiento para el hallazgo de la solucion del problema 0OJ - HH aqui estudiado, incluye el 
caso de su version no singular. En efecto, sin perder generalidad, supongamos que en la ecuacion 

f K 



N 



J r -1 , la matriz 11 + A 



T 1 . Escogiendo A suficientemente grande, esto 



(1401 , E = I. Al considerar la descomposicion canonica de A = T 

1T , T / XI + K 
obtiene la forma Al + A = T n , T AT 

\ A.I + N 

hara que la diagonal de los bloques XI + K y A7 + N no se anule y asi, la condition: que existe y e C 
tal que {yE + A) es no singular, es verificada, y nuestro procedimiento puede ser aplicado. Notese 
ademas que E E = Iy E A = A, por lo tanto la condition 1T75I I puede omitirse. 
2. Si la condition 05J del teorema |22] no se cumple, entonces de acuerdo a la proposicion |20j la 
ecuacion J59j puede ser reducida a lEHJ ; como la solucion de esta es la solucion trivial, se concluye 
que asi tambien lo sera la solucion encontrada para J50j - J54J . 



14 



MANUEL J. SALAZAR, EDISON E. VILLA 



4.4. Conclusion. En este trabajo se ha demostrado la existencia de soluciones numericas estables pa- 
ra una gama de sistemas de ecuaciones en derivadas parciales singulares de la ecuacion de onda ho- 
mogenea como la descrita en lO - dm. En nuestra metodologia, con respecto a los eigenvalores del 
problema de Sturm Liouville discreto asociados a Q3J - CES), hemos asumido el teorema espectral para 
matrices simetricas en vez de aplicar metodos numericos para establecer su existencia. Dicho enfoque 
nos permitio enfrentar con otra estrategia el problema de encontrar tal solution a nuestro sistema de 
interes para el cual la arbitrariedad de los parametros a y ft, relacionados al PSD {T4J - {T6J puede ser 
usada con el fin de que las hipotesis del teorema (122) se verifiquen. 
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